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Abstract 


Studies of rest-frame optical emission in quasars at z>6 have historically been limited by the wavelengths 
accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to 
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probe this emission deep into the reionization epoch. We report the observations of eight quasars at z > 6.5 using 
the JWST/NIRCam Wide Field Slitless Spectroscopy as a part of the “A SPectroscopic survey of biased halos In 
the Reionization Era (ASPIRE)" program. Our JWST spectra cover the quasars' emission between rest frame 
~4100 and 5100 A. The profiles of these quasars’ broad H8 emission lines span a full width at half maximum from 
3000 to 6000 km s^ !. The H-based virial black hole (BH) masses, ranging from 0.6 to 2.1 billion solar masses, 
are generally consistent with their Mg II-based BH masses. The new measurements based on the more reliable HG 
tracer thus confirm the existence of a billion solar-mass BHs in the reionization epoch. In the observed [O m] AA 
4960,5008 doublets of these luminous quasars, broad components are more common than narrow core components 
(€ 1200 km s^), and only one quasar shows stronger narrow components than broad. Two quasars exhibit 
significantly broad and blueshifted [O m] emission, thought to trace galactic-scale outflows, with median velocities 
of —610 and —1430 km s! relative to the [C II] 158 jum line. All eight quasars show strong optical Fe II emission 
and follow the eigenvector 1 relations defined by low-redshift quasars. The entire ASPIRE program will eventually 
cover 25 quasars and provide a statistical sample for the studies of the BHs and quasar spectral properties. 


Unified Astronomy Thesaurus concepts: Quasars (1319); Supermassive black holes (1663); Reionization (1383) 


1. Introduction 


Observations of luminous z = 6.5 quasars have revealed the 
existence of supermassive black holes (SMBHs) and their 
massive host galaxies in the reionization epoch. The measure- 
ments of high-redshift quasar BH masses suggest that 105? 
solar-mass SMBHs existed as early as 700 million yr after the 
Big Bang, challenging models of early SMBH growth and BH 
seed formation (e.g., Bafiados et al. 2018; Yang et al. 2020; 
Wang et al. 2021; Yang et al. 2021; Farina et al. 2022). 
Measuring BH masses accurately is the key step required to 
characterize the growth and evolution of early SMBHs. The 
way to obtain robust BH mass measurements is reverberation 
mapping (RM). However, this is not practically possible for 
z>4 quasars. At higher redshift, the BH mass estimates are 
based on the so-called single-epoch virial BH mass using 
quasar broad emission lines. 

The observations of quasar spectra in the rest-frame optical 
have been mostly confined to quasars at z < 4 (e.g., Shen 2016; 
Matthews et al. 2021) because of the difficulty in observing 
faint objects from the ground at wavelengths longer than the 
near-infrared K band. Thus, the measurements of SMBHs in 
quasars at z > 4 are mainly based on broad emission lines in the 
rest-frame UV (e.g., CIV A1549 and Mg II \2800). The Mg II- 
based BH mass estimators are widely used for the measure- 
ments of reionization-era SMBHs (e.g., Shen et al. 2019; Yang 
et al. 2021; Farina et al. 2022). However, the study of quasar 
BH masses at low redshift via RM and the comparison between 
multiple single-epoch BH mass estimators have established the 
HQ line as the most reliable tracer among quasar UV /optical 
emission lines (e.g., Shen 2013; Wang et al. 2020). The 
existing Mg Il-based scaling relation is indirectly calibrated 
based on the H relation derived from RM, which introduces 
additional uncertainties. Therefore, high-sensitivity infrared 
(IR) spectroscopic observations covering the rest-frame optical 
emission are needed to better estimate the HG-based BH masses 
of quasars at the highest redshift. 

A SPectroscopic survey of biased halos In the Reionization 
Era (ASPIRE) is a James Webb Space Telescope (JWST) 
Cycle 1 GO program that will observe a sample of 25 z > 6.5 
quasars with JWST/NIRCam imaging and Wide Field Slitless 
Spectroscopy (WESS; Wang et al. 2023) The WFSS 
observations will provide the spectra in the H8 region for 
reionization-era quasars for the first time. In addition, the high- 
quality JWST spectra will allow us to study the spectral 
properties (e.g., continuum emission; [O m] AA 4960, 5008; 
and optical Fe II emission) of quasars in the early universe. In 


this paper, we report JWST observations of eight ASPIRE 
quasars. They form the first sample of reionization-era quasars 
with spectral coverage in the rest-frame optical. We describe 
the observations and data reduction in Section 2. In Section 3, 
we report the main results from this early data, including the 
spectral fitting procedure, the HG-based BH masses, and the 
quasar optical emission-line properties. In Section 3, we also 
present detailed analysis of the blueshifted [O Mm] components 
and outflows observed in these quasars. We summarize our 
work in Section 4. All results assume a flat ACDM cosmology 
with parameters Qa = 0.7, Qm = 0.3, and hA — 0.7. 


2. Observations and Data Reduction 


The eight quasar spectra used in this work are part of the 
medium-size JWST Cycle 1 program ASPIRE (PID 2078, PI: 
F. Wang; Wang et al. 2023), which carries out a spectroscopic 
survey in 25 reionization-era quasar fields using JWST/ 
NIRCam WFSS using the F356W filter, as well as imaging 
observations in the F115W, F200W, and F356W bands at the 
same time. This program allows us to statistically investigate 
both reionization-era quasars and high-redshift galaxies in the 
quasar fields. The NIRCam spectral and imaging data set 
probes the rest-frame optical emission in a large sample of 
z> 6.5 quasars and host galaxies for the first time, which will 
improve our understanding of early SMBH activity and their 
host galaxies. The spectroscopic and imaging observations of 
the quasar fields allow studies of the environments of 
reionization-era quasars in a statistical manner, which will 
address the long-standing question of whether luminous 
quasars reside in the most massive dark matter halos and 
inhabit large-scale galaxy overdensities in the early universe. 

The eight quasars reported in this work were observed with 
JWST/NIRCam in 2022 August and September. These targets 
were observed using the R-grism with F356W filter in long 
wavelength (LW) and F200W imaging in short wavelength 
(SW) simultaneously. We performed a three-point INTRAMO- 
DULEX primary dither with a 2-POINT-LARGE-WITH- 
NIRISS subpixel dither. A SHALLOW4 readout mode was 
used with nine groups per integration and one integration per 
exposure, which resulted in a total exposure time of 2834 s. 
After the grism observation, direct image and out-of-field 
images were taken in F356W (LW) and F115W (SW) with a 
total exposure time of 1417 s. For the analysis in this work, we 
only use the WESS F356W spectral data of the quasar targets. 
Analysis of the spectra of galaxies in the quasar fields and 
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broadband imaging data will be presented in a set of 
subsequent papers. 

The data reduction of imaging and WFSS data will be 
presented by Wang et al. (2023) in detail, and we will briefly 
describe the related steps below. The reduction of WFSS data 
used in this work was based on the JWST pipeline 1.7.2 and 
calibration reference data “jwst_0989.pmap.” We started with 
calwebb detectorl to process the raw data and applied a 
custom step to * rate.fits to remove the 1/f noise feature 
(Schlawin et al. 2020). We performed flat-fielding and assigned 
World Coordinate System for exposures and then generated 
master background models based on all existing ASPIRE 
observations. We constructed a tracing model and a dispersion 
model following the method used in Sun et al. (2022) and 
applied 2D spectral extraction to background-subtracted data. 
We used updated sensitivity functions based on Cycle 1 
calibration programs (PID 1536, 1537, and 1538). Then the 1D 
spectra were extracted from individual exposures using optimal 
extraction (Horne 1986) and coadded with inverse variance 
weighting. The 1D spectra of the eight ASPIRE quasars are 
shown in Figure 1. 


3. Rest-frame Optical Spectra of Reionization-era Quasars 


In this section, we describe our spectral fitting procedure and 
the main results derived from this early sample. We measure 
the BH masses based on the H8 line width and continuum 
luminosity and compare them with the measurements derived 
from their rest-frame UV spectra. We also discuss the profile of 
the [O mM] emission lines and search for potential galactic-scale 
outflows. We note that the data reduction procedure of this new 
instrument is still in development, and the spectra might be 
slightly affected by future changes of the reduction pipeline, 
but the main findings in this paper will not be significantly 
affected. The final JWST ASPIRE data set will be presented 
after this program has been completed. 


3.1. Spectra Fitting 


We perform spectral fitting for each WFSS spectrum using a 
model consisting of continuum and emission lines. The redshift 
derived from the [C II] 158 uum line (Venemans et al. 2020; 
Yang et al. 2021; Wang et al. 2023) is used as the initial 
redshift. We use spectral data between observed wavelengths 
31600 and 39500 A, excluding edge regions where the 
calibration is less certain. This covers the rest-frame wave- 
length range of —4100—5200 A for these quasars. We use a 
pseudocontinuum model including a power-law continuum and 
an empirical optical Fell template from Boroson & Green 
(1992). To fit this pseudocontinuum, we choose continuum 
windows that are free of strong emission lines in quasar 
composite, [4150, 4230], [4435, 4700], and [5100, 5200] A in 
the rest frame. For the two highest-redshift quasars, 
J0109—3047 and J0218--0007, we adjust the windows due to 
their limited spectral coverage on the red side. For the former, 
we only use the first two windows. For the latter, we extend the 
blue-side cutoff of the last window to 5060 A. 

H6 and [Om] lines. We adopt the rest-frame vacuum 
wavelengths of 4862.68 A for the Hø line and 4960.30 A and 
5008.24 A for the [O m] doublet (Vanden Berk et al. 2001). We 
employ a line model consisting of one narrow Gaussian 
component and three broad Gaussian components for the H8 
line (e.g., Shen 2016). For each [Om] line, we use one 
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Gaussian for the narrow component (the “core” component) 
plus a second Gaussian profile for the potentially broadened 
and blueshifted component. For the object with strong 
blueshifted [O m] wings (i.e., J0224-4711), we include one 
more Gaussian for its broad [O III] component. We set a lower 
limit to the line FWHM of each component to be 200 km s~". 
This is based on the instrument resolution (R ~ 1600 at 4 um). 
We also set an upper limit of 1200 kms ' for all narrow 
components (e.g., Shen et al. 2011). 

We tie the velocity shift and dispersion of the narrow H to 
that of the [O m] “core” component. We assume a fixed line 
ratio of 3.0 between [O m] A5008 A and [O m] A4960 A for 
both "core" and broad components (e.g. Osterbrock & 
Ferland 2006). The velocity shifts of the broad [O m] 
components of each transition are also tied. We do not fix 
the line ratio between narrow H8 component and narrow [O III] 
lines but limit the ratio of the “core” [O m] 5008 to the narrow 
HÓ to Z1. If for an object, there is no narrow [O m] 5008 line, 
we consider that there should also be no narrow Hf. For 
objects with narrow [O III] components, we derive [O III]-based 
redshifts using their narrow components and find that the 
[O m]-based redshifts show good agreement with their redshifts 
derived from [C II] 158 um lines. 

He line. We test the possibility of including a He II line 
with a vacuum wavelength of 4687.02 A by running alternate 
spectral decomposition. We include one Gaussian profile 
without constraining the line width and use a continuum 
window of [4435, 4635] A instead of [4435, 4700] A to avoid 
potential contamination from He II emission. We run the same 
spectral fitting procedure and fit the He II line together with the 
H@ and [Om] lines. We find that only one quasar (J2232 
+2930) has an obvious Hell line. Another quasar (J0218 
+0007) might have a narrow Hell line, but it only has a 
significance of ~2ø. For these two, including He II reduces 
their BH masses by 4% and 2%, respectively. For other objects, 
their best-fit models do not show any need for a Hell line, 
while the shorter continuum windows affect fits of their 
continuum, especially for z > 6.7 quasars. Thus, we finally only 
apply the He II fit for quasar J2232--2930, and for other quasars 
we keep the continuum window cutoff of 4700 A. Even if there 
could be some weak He II emission in these quasars, their BH 
masses would not be affected significantly. 

Hy line. We also apply a two-Gaussian fit to the Hy line at 
4341.68 A. The two Gaussian profiles are set for one narrow 
and one broad component. The same limits of narrow line 
width as used for H are applied. The results are only used for 
plotting the best-fit but not any science analysis in this work. 

We utilize a Monte Carlo (MC) approach to estimate the 
uncertainties of all spectral measurements, following Yang et al. 
(2021; also Shen et al. 2019; Wang et al. 2020). For each 
spectrum, we generate 50 mock spectra by randomly adding 
Gaussian noise at each pixel with a standard deviation equal to the 
spectral error at that pixel. We then perform the same spectral 
fiting procedure to each mock spectrum. The uncertainty of each 
spectral measurement is estimated by averaging the 16% and 84% 
percentile deviations from the median value. The best-fit results of 
all quasars are shown in Figure 1. 


3.2. HB-based Black Hole Masses 


The single-epoch virial method is widely used for estimating 
the BH mass of high-redshift quasars, assuming virial motion of 
the line-emitting gas in the quasar broad-line region (BLR) and 
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Figure 1. The JWST NIRCam WFSS F356W spectra of eight ASPIRE quasars (black line) with spectral uncertainty (gray), ordered by R.A. The red solid lines denote 
the best total fits. The best fits of different spectral components are shown with dashed and dotted lines. For each object, the left panel presents the entire spectrum with 
total fit, power-law continuum, and pseudocontinuum (i.e., power-law plus Fe II emission). The right panel shows the zoomed-in region of H8 and [O Il] lines. The 
systemic redshifts of the quasars are based on their [C II] line measurements from ALMA observations. The blue squares with error bars are the photometric data in the 
WISE W1 3.4 um band. The panels below the spectra present the residuals (data model) of each best-fit model. 


applying the empirical correlation between the BLR size and 
quasar continuum luminosity (i.e., the A—L relation). With this 
method, we can estimate the quasar central BH mass utilizing 
quasar broad emission lines and continuum luminosity. Quasar 
UV and optical emission lines, C IV A1549 A, Mg I 42800 A, Ha 


6563 A, and Hf have all been used as virial BH mass tracers 
(e.g., McLure & Jarvis 2002; Vestergaard 2002; McLure & 
Dunlop 2004; Onken et al. 2004; Greene & Ho 2005; Vestergaard 
& Peterson 2006; Vestergaard & Osmer 2009; Shen et al. 2011). 
Calibration coefficients used in these BH mass estimators are 
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Figure 1. (Continued.) 


determined using low-redshift quasar samples that have BH mass 
measurements based on RM. Among these lines, H8 is often 
thought to be the most reliable tracer from both RM-based 
measurements and comparison of single-epoch mass estimators 
(e.g., Shen & Liu 2012). 

We measure the single-epoch virial BH masses of these 
quasars based on their continuum luminosity at rest-frame 
5100 A and the H@ line width, obtained from the spectral 


fitting. We adopt the BH mass estimator in Vestergaard & 
Peterson (2006), 
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Msau,H l 
—s ad (oe 


Mo 1000 km s^! 10% erg s~ 


The 5100 À continuum luminosity is derived from the best-fit 
power-law continuum directly. The systematic scatter of this 
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Figure 1. (Continued.) 


BH mass estimator relative to RM mass is ~0.4—0.5 dex, which 
is not included when we report the errors of BH mass 
measurements. The uncertainties of all measurements are 
derived from the MC approach, as described above. We then 
compute the Eddington ratio (Agaa = Lpo1/Leaa) for each 
quasar. The bolometric luminosity is estimated by multiplying 
the 5100 A luminosity by 9.26 (Richards et al. 2006). Table 1 
summarizes the redshift, H@ line width, continuum luminosity, 
HG BH mass, Eddington ratio, and Mg II BH mass of each of 
the eight quasars. 

The H6-based measurements show that these quasars have 
BH masses in the range of (0.6-2.1) x 10? Mo. These quasars 
all have been published with measurements of BH mass based 
on the MgII emission lines from their near-infrared (NIR) 
spectroscopic data (e.g., Mazzucchelli et al. 2017; Reed et al. 
2019; Schindler et al. 2020; Wang et al. 2021; Yang et al. 
2021; Farina et al. 2022). We adopt the most recent 
measurements of these eight quasars from Reed et al. (2019), 
Wang et al. (2021), Yang et al. (2021), and Farina et al. (2022). 
All of these Mg II-based BH masses are computed using the 
BH mass estimator from Vestergaard & Osmer (2009). 
Figure 2 shows the Hó-based BH masses compared to 
their MgII-based BH masses. The BH masses derived from 
the H@ lines and 5100 A luminosity are consistent with the 


Mg II-based measurements within their systematic scatters (i.e., 
7-0.4—0.5 dex). Recently, a new measurement of BH mass of 
quasar J01004-2802 at z = 6.3 using the HØ line from JWST 
observations also shows agreement with its Mg II-based BH 
mass within systematic uncertainties (Eilers et al. 2022). 

In our quasar sample, the differences between these two BH 
mass estimates span from —0.32 to 0.22 dex, comparable to the 
scatter shown in the SDSS low-redshift quasar sample (a 
standard deviation of —0.38 dex; Shen et al. 2011). Among 
these eight quasars, there is a trend that the H6-based BH 
masses are slightly smaller, with an average Mpu,uo-MnpnMg u 
of —0.13 dex. However, the current sample is too small to 
conclude if there is a systematic offset. The H-based BH 
masses and the comparison between H8 and Mg II measure- 
ments confirm the existence of a billion solar-mass BHs in the 
reionization epoch. These SMBHs require za Mo seed 
BHs at z — 30 assuming Eddington accretion across the entire 
time since BH seeding and a radiative efficiency of 0.1 (e.g., 
Yang et al. 2021; Farina et al. 2022). This thus raises the 
questions of how to maintain highly accreting BHs for the long 
term and how to form massive seed BHs in the early universe 
(e.g., Volonteri 2012; Davies et al. 2019; Inayoshi et al. 2020). 

These measurements are derived from the spectral fits of the 
NIR spectra and JWST spectra separately. Thus, for both data 
sets, the limited continuum windows might lead to additional 


Table 1 

Spectral Fitting and Quasar Properties 
Name Z[c m] Zfo m? Manus FWHMmg; Ls100 Azad Mpue 1 Ref Mg 1” 

(0? Mo) (km s^!) (109erg s^!) (10° Mo) 
J010953.13—304726.30 6.7904 + 0.0003 e 0.60 + 0.05 3033.0 + 117.0 6.40 + 0.08 0.76 + 0.06 1.11704, Farina+2022 
J021847.04--000715.20 6.7700 + 0.0013 6.7668 + 0.0006 0.61 + 0.05 3030.0 + 117.0 6.71 + 0.16 0.78 + 0.07 0.61 + 0.07 Yang+2021 
J022426.54-471 129.40 6.5222 + 0.0001 6.5186 + 0.0004 2.15 + 0.09 3936.0 + 78.0 29.19 + 0.09 0.97 + 0.04 1.30 + 0.18 Wang+2021 
J022601.874-030259.28 6.5405 + 0.0001 e 1.05 0.09 3131.0 + 129.0 17.25 + 0.09 1.17 + 0.09 2.20 + 0.39 Wang+2021 
J024401.02—500853.70 6.7306 + 0.0002 6.720 + 0.003 0.88 + 0.10 2969.0 + 158.0 15.08 + 0.13 1.22 + 0.12 1.15 + 0.39 Reed+2019 
J030516.92-315056.00 6.6139 + 0.0002 6.615 + 0.001 0.62 + 0.14 3020.0 + 314.0 6.95 + 0.06 0.80 + 0.14 0.70 + 0.38 Wang+2021 
J200241.59-301321.69 6.6876 + 0.0004 Us 0.84 + 0.04 2981.0 + 79.0 13.65 + 0.10 1.15 + 0.06 1.62 + 0.27 Yang+2021 
J223255.15+293032.04 6.666 + 0.004 6.669 + 0.0001 1.93 + 0.13 6044.0 + 201.0 4.24 + 0.10 0.16 + 0.01 3.06 + 0.36 Yang+2021 
Notes. 


* The [O n1] redshift is estimated based on the line center of the “core” component of [O m] 5008 line. Quasars J0109—3047, J0226--0302, and J2002-3013 do not have narrow [O III] components, and thus we do not 
estimate their [O III] redshift. 
> The references of the Mg II-based BH masses adopted in this work. 
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Figure 2. The BH masses of eight ASPIRE quasars (red squares), measured 
from the H£ lines in JWST/NIRCam WFSS spectra and the Mg II lines in their 
NIR spectra. These are compared with the measurements of low-z quasars from 
the SDSS Data Release 7 (DR7) sample (Shen et al. 2011) using the same 
Mg II and H8 BH mass estimators (gray small squares and blue contours). The 
differences between Mgy,ng and Mgyme n Of these high-redshift quasars are 
comparable to the scatter shown in the low-redshift quasar sample. The mean 
value of Mau ua-Mpn Mg n is —0.13 dex, while the sample size is too small to 
identify a systematic offset. 


uncertainties of the spectral decomposition and BH mass 
measurements. In the next step, we will utilize the entire 
ASPIRE sample to perform a detailed comparison between the 
Hf and Mg II BH masses when observations of the full sample 
are complete. With the final sample, we will combine the NIR 
and JWST spectra as well as the JWST F115W, F200W, and 
F356W photometries and will apply joint spectral analysis. The 
combined data set will provide a wider continuum window and 
allow a better separation of the Fell emission from the 
featureless continuum. 


3.3. Optical Spectral Properties and [O iti] Outflows 


The most important spectral quantities covered by the WFSS 
spectra are the H@, [O rri], and Fe II emission. In general, all 
eight quasars have strong and broad Hj lines, with FWHM 
ranging from 3000 to 6000 km s '. Optical Fe II emission is 
clearly seen in all eight quasars. The observed [O mM] doublets 
in these quasars have relatively diverse properties. Three of 
these quasars, J0109—3047, J02264-0302, and J2002-3013, 
show similar profiles in their [O III] line regions, with very 
weak broad-only [O mM] doublets but no narrow [O m] lines ( 
ie. the "core" component, «1200 km s^). The other five 
quasars all have both narrow and broad [O m] components. 
Quasar J2232+2930 has the strongest narrow [Om] lines 
among these quasars. It is the only object observed with 
stronger narrow [O IN] lines relative to its broad components. 

Using the WFSS data, we are able to visit the eigenvector 1 
(EV1) relations (Boroson & Green 1992) in the most distant 
quasars. The EV1 relations describe the relationships between 
emission-line quantities and the strength of the optical Fe II 
emission and has been widely investigated in the past three 
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Figure 3. EV1 relations using our eight high-redshift quasars (red squares), 
compared with luminous quasars at intermediate redshift (Shen 2016; open 
diamonds) and SDSS low-z quasars (Shen et al. 2011; gray dots). The black 
filled circles and the dashed lines denote the median, 16th, and 84th percentiles 
of the SDSS low-z quasars at each redshift bin. These reionization-era quasars 
generally follow the trends observed in lower-redshift quasars. 


decades with quasar samples at z«4 (e.g, Boroson & 
Green 1992; Marziani et al. 2001; Boroson 2002; Shen & 
Ho 2014). The anticorrelation between the strength of [O m] 
and Fe II emission is one of the most important EV1 relations 
observed. The trend of Fell strength, as the HG FWHM 
decreases, is also a probe of the physical drivers of the EV1 
relations. We study the correlation of the [O m] 5008 EW and 
FWHMy, with the Rg, n (Shen 2016) and compare them with 
the observations of the luminous intermediate-redshift 
(1.5 <z<3.5) sample in Shen (2016) and the low-redshift 
(z « 1) quasars from the SDSS DR7 (Shen et al. 2011). The 
EWio m5008 used here is the EW of the entire [O m] A5008 line 
in order to compare our results with the measurements from 
Shen (2016). FWHMgj; is the line width of the broad HØ. Rp, n 
is the ratio of EWzg, u and EWy;, and the Fe II EW is measured 
between 4434 and 4684 A. 

As shown in Figure 3, overall, the optical spectral quantities 
in these high-redshift quasars follow the relations defined in 
lower-redshift quasars. Our high-redshift quasars show slightly 
lower [O m] EW than the mean values of SDSS low-redshift 
quasars and are more similar to the luminous 1.5 «z «3.5 
quasar sample, which is due to their high luminosity and the 
Baldwin effect (Baldwin 1977). There is no significant outlier 
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among our high-redshift quasars. The three lowest [O m] EWs 
are from the quasars with weak and broad-only [O II] emission. 
In the FWHMygg versus Ry, n plane, our objects distribute in a 
narrow FWHMy, range with lower FWHMy, than the 
intermediate-redshift sample. This is consistent with the high 
Eddington ratios of most (7 of 8) of our high-z quasars. The 
general consistency between our sample and low-z distributions 
suggests that the EV1 relations might exist in quasars as early 
as the reionization epoch, and our full ASPIRE sample will 
enable more detailed comparisons with low-z samples. 

Broad, skewed, and blue/redshifted [O m] profiles have 
been discovered in quasars at z<3.5 and used to probe 
kiloparsec-scale outflows in quasar host galaxies (e.g., Cano- 
Diaz et al. 2012; Carniani et al. 2015; Zakamska et al. 2016; 
Bischetti et al. 2017). At z > 6, outflows powered by the central 
active galactic nucleus are a common prediction of a broad 
range of cosmological simulations (e.g., Dubois et al. 2013; 
Costa et al. 2014, 2022; Lupi et al. 2022). However, 
observational evidence of large-scale outflows in quasar hosts 
remains anecdotal at this redshift. The JWST data set provides 
us with the first chance to search for [Om] outflows in 
reionization-era quasars. Two quasars, J02184-0007 and 
J0224—4711, have prominently strong broad and blueshifted 
[O m] components, with a velocity shift of —630 kms! and 
—1690kms ! relative to the narrow [Om], respectively. 
Figure 4 shows their zoomed-in [O mM] 45008 lines, with other 
components subtracted out according to the best-fit models. We 
compute the median velocity, vso, the velocity at the 50th 
percentile of the line flux; vio and vog, the velocities at the 10th 
and 90th percentiles; and wso (—voo — v10), containing 80% of 
the total line power, which is close to the commonly used 
FWHM (wgo— 1.088 FWHM) for a Gaussian profile. The 
velocities are calculated relative to their [C II] redshifts. J0218 
+0007 has a median velocity (vso) of ~—610 kms, with a 
voo of —1510 km s^! and a wep of 1590 km s~". 0224-4711 
has a median velocity (vso) of ——1430 km s ' a vo of 
—4010 km s™!, and a wgọ of ~4140 km s™!. 

Their broad [O mM] widths and high-velocity shifts make 
them stand out among the quasar population. Particularly, 
J0224—4711 has one of the most extreme broad and blueshifted 
[O m] lines observed to date, even compared to the observa- 
tions of lower-redshift quasars. For example, the 1.5 < z < 3.5 
luminous quasars in Shen (2016) have a median [O m] FWHM 
of 900kms !, with the broadest one to 2100 km s^ !. 
Zakamska et al. (2016) observed extreme outflows in red 
luminous quasars at z ~ 2-3 and discovered broad [O III] with 
Wgo ranging between 3600 and 5500 kms !, strongly blue- 
shifted up to 1500km s |. The more recent study of 
hyperluminous quasars at z~ 2-4 revealed a population of 
[O m] lines with FWHM of 1200-2200 km gn (Bischetti et al. 
2017). Large velocity widths and high velocities of blueshifted 
[O mM] suggest powerful ionized outflows, while the nature of 
high-velocity large-scale gas in quasars is still debated. The 
most blueshifted [O mM] lines observed in the extremely red 
z = 2-3 quasars are thought to be caught during a “blow-out” 
phase of quasar evolution and associated with fast ionized 
outflows (e.g., Zakamska et al. 2016; Perrotta et al. 2019; 
Vayner et al. 2021). So far, there is no evidence of J0218 
+0007 and JO224-4711 being extremely red quasars. A 
detailed multiwavelength analysis will be conducted in the 
future to better understand the physical properties of their fast 
outflows. In addition, a joint analysis of the 1D spectra and 2D 
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Figure 4. The [O m] 5008 lines with strong broad and blueshifted components 
have been observed in two quasars, J0218--0007 (top) and J0224—-4711 
(bottom). The blue solid lines denote the median velocity (vso), and the dashed 
lines mark the part containing 80% of the line power (wgo). J02184-0007 has a 
vso of ~—610 km s^! and a wgo of 1590 km s^ '; quasar J0224—4711 has a 
more remarkable [O m], with a vsọ of ~—1430kms ! and a weg of 
4140 km s^!. These blueshifted broad lines, in particular in J0224—4711, are 
comparable to the most extreme cases observed at lower redshift (e.g., 
Zakamska et al. 2016; Bischetti et al. 2017). 


images will be performed in subsequent works to examine the 
host galaxies and outflows. In the next step, JWST NIRSpec/ 
IFU follow-up observations will be carried out to reveal the 
nature of their outflows. The entire ASPIRE sample will allow 
us to systematically study [O m] outflows in the reionization- 
era quasars, in synergy with multiwavelength imaging and 
spectroscopic data. 


4. Summary 


In this paper, we present the results from an early sample of 
ASPIRE JWST/WFSS spectral data set consisting of eight 
quasars at z> 6.5, observed in 2022 August and September. 
With these high-quality WESS spectra, we are able to study the 
rest-frame optical emission in reionization-era quasars for the 
first time. We measure the virial BH masses using the H£ lines 
and compare them with the Mg II-based measurements. We 
find that, in general, the HG BH masses are consistent with the 
Mg II BH masses within the systematic uncertainties of the BH 
mass estimators. The discrepancy is within the scatter between 
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these two measurements as gauged from lower-redshift quasar 
samples. Therefore, using a more reliable tracer, the H5 line, 
we confirm the existence of a billion solar-mass BHs in the 
reionization epoch. 

We also test the optical spectral properties of these 
reionization-era quasars in the EV1 plane defined by z «4 
quasars. The results suggest that EVI relations exist in these 
distant quasars and our quasars are distributed similarly to the 
luminous quasars at intermediate redshift. We discover two 
quasars with broad and strongly blueshifted [O m] lines. One of 
them has a prominent broad [O M] profile (wg9 = 4140 km s!) 
with a high-velocity shift (vso = —1430 km s |, which is 
among the most extreme cases, even compared to observations 
in lower-redshift quasars. 

With the ASPIRE program, we will construct a sample of 25 
quasars with not only JWST WFSS spectra and NIRCam 
images but also the data from multiwavelength observations, 
spanning X-ray to submillimeter and radio. The final sample 
will be used for detailed investigations of their BH physics, 
quasar feedback, and host galaxy properties. 
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